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Brane world cosmologies seem to provide an alternative explanation for the present accelerated 
stage of the Universe with no need to invoke either a cosmological constant or an exotic quintessence 
component. In this paper we investigate statistical properties of gravitational lenses for some par- 
ticular scenarios based on this large scale modification of gravity. We show that a large class of 
such models are compatible with the current lensing data for values of the matter density parameter 
£3 m < 0.94 (ler). If one fixes fi m to be ~ 0.3, as suggested by most of the dynamical estimates of 
the quantity of matter in the Universe, the predicted number of lensed quasars requires a slightly 
open universe with a crossover distance between the 4 and 5-dimensional gravities of the order of 
1.76-H7 1 . 

PACS numbers: 98.80.Es; 04.50.+h; 95.35.+d 



I. INTRODUCTION 



The results of observational cosmology in the last years 
have opened up an unprecedented opportunity to test the 
veracity of a number of cosmological scenarios as well 
as to establish a more solid connection between parti- 
cle physics and cosmology. The most remarkable finding 
among these results comes from distance measurements 
of type la supernovae (SNe la) that suggest that the ex- 
pansion of the Universe is speeding up, not slowing down 
jjj . As widely known such a result poses a crucial prob- 
lem for all CDM models since their generic prediction 
is a decelerating universe (q Q > 0), whatever the sign 
adopted for the curvature parameter. Indirectly, similar 
results have also been obtained, independent of the SNe 
la analyses, by combining the latest galaxy clustering 
data with CMB measurements 0] . 

To reconcile these observational results with theory, 
cosmologists have proposed more general models con- 
taining a negative-pressure dark component that would 
be responsible for the present accelerated stage of the 
Universe. Although a large number of pieces of obser- 
vational evidence have consistently suggested a universe 
composed of ~ 2/3 of dark energy, the exact nature of 
this new component is not well understood at present. 
Among the several candidates for dark energy discussed 
in the recent literature, the simplest and most theoreti- 
cally appealing possibility is the vacuum energy or cos- 
mological constant. Despite the serious problem that 
arises when one considers a nonzero vacuum energy 
models with a relic cosmological constant (ACDM) seem 
to be our best description of the observed universe, being 
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considered as a serious candidate for standard cosmology. 

On the other hand, motivated by particle physics con- 
siderations, there has been growing interest in cosmo- 
logical models based on the framework of brane-induced 
gravity j|, ||, 0|. The general principle behind such 
models is that our 4-dimensional Universe would be a 
surface or a brane embedded into a higher dimensional 
bulk space-time on which gravity can propagate. In some 
of these scenarios, there is a certain crossover scale r c that 
defines what kind of gravity an observer on the brane will 
observe. For distances shorter than r c , such an observer 
will measure the usual 4-dimensional gravitational 1 /r 2 
force whereas for distances larger than r c the gravita- 
tional force follows the 5-dimensional 1/r 3 behavior. In 
this way, gravity gets weaker at cosmic distances and, 
therefore, it is natural to think that such an effect has 
some implications on the dynamics of the Universe S . 

Several aspects of brane world cosmologies have been 
explored in the recent literature. For example, the issue 
related to the cosmological constant problem has been 
addressed as well as evolution of cosmological per- 
turbations in the gauge- invariant formalism flCfl , cosmo- 
logical phase transitions pl[ , inflationary solutions [ OJ , 
baryogenesis Jll| , stochastic background of gravitational 
waves pi , singularity, homogeneity, flatness and entropy 
problems [T|] , among others (see [|16[ for a discussion on 
the different perspectives of brane world models) . From 
the observational viewpoint, however, the present situ- 
ation is somewhat controversial. While the authors of 
Refs. [[ijj |l8) have shown that such models are in agree- 
ment with the most recent cosmological observations (for 
example, they found that a flat universe with Q m = 0.3 
and r c ~ 1AH~ is consistent with the currently SNe la 
+ CMB data) , the authors of Ref. |ll| have claimed that 
a larger sample of SNe la data can also be used to rule 
out these models at least at the 2a level. Recently, one of 
us 0] used measurements of the angular size of high-z 
compact radio sources to show that the best fit model for 
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these data is a slightly closed universe with f2 m ~ 0.06 
and a crossover radius of the order of 0.94_ff o _1 . 

For the reasons presented earlier, the comparison be- 
tween any alternative cosmology and ACDM models is 
very important. In this concern, statistical properties of 
gravitational lenses may be an interesting tool because, 
as is well known, they provide restrictive limits on the 
vaccum energy contribution (see, for instance, [pl|). On 
the other hand, in brane world models the distance to an 
object at a given redshift z is smaller than the distance 
to the same object in ACDM models (assuming the same 
value of f2 m ). Therefore, we expect that the constraints 
coming from lensing statistics will be weaker for these 
models than for their ACDM counterparts. 

In this paper, we explore the implications of gravita- 
tionally lensed QSOs for models based on the framework 
of the brane-induced gravity of Dvali et al. || that have 
been recently proposed in Refs. || [|. We restrict our 
analysis to accelerated models, or equivalently, models 
that have a "self- inflationary" solution with H ~ r^ 1 (H 
is the Hubble parameter). As explained in jfj), in such 
scenarios, the bulk gravity sees its own curvature term 
on the brane as a negative-pressure dark component and 
accelerates the Universe. 

This paper is organized as follows. In section 2, we 
present the basic field equations and distance formulas 
relevant for our analysis. We then proceed to analyze 
the constraints from lensing statistics on these models 
in section 3. In section 4 our main conclusions are pre- 
sented. 



II. THE MODEL: BASIC EQUATIONS AND 
DISTANCE FORMULAS 



The Friedmann's equation for the kind of models we 
are considering is |8[ p7[ 
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where p is the energy density of the cosmic fluid, k = 
0, ±1 is the spatial curvature, M p i is the Planck mass 
and r c — M^/2M§ is the crossover scale defining the 
gravitational interaction among particles located on the 
brane (M5 is the 5-dimensional reduced Planck mass). 
From the above equation we find that the normalization 
condition is given by 



y/n Ic + n m = 1 



(2) 



where Q m and fl^ are, respectively, the matter and cur- 
vature density parameters (defined in the usual way) and 



l/4r 2 c H 2 , 



(3) 



becomes ITi 



for f2 rc < 1 and fi ro < 1. (4) 



In order to derive the constraints from lensing statistics 
in the next section we shall use the concept of angular 
diameter distance, Da(z). Such a quantity can be eas- 
ily obtained in the following way: consider that photons 
are emitted by a source with coordinate r = n at time 
ti and are received at time t a by an observer located at 
coordinate r = 0. The emitted radiation will follow null 
geodesies on which the dimensionless comoving coordi- 
nates 6 and <j> are constant. The comoving distance of 
the source is defined by 



dt 



dR 
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From Eqs. (1) and (5), it is possible to show that the 
comoving distance r\(z) can be written as p0| 
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where the subscript o denotes present day quantities, 
x' = Ijp- = (l + z) -1 is a convenient integration variable 
and the function Yl ( r ) ls defined by one of the follow- 
ing forms: XX r ) = sinh(r), r, and sin(r), respectively, 
for open, flat and closed geometries. The dimensionless 
function f{Slj, x) is given by 



1/2 



(7) 



where j stands for m, r c and k. 

The angular diameter distance to a light source at 
r = r\ and t = t\ and observed at r = and t = t 
is defined as the ratio of the source diameter to its angu- 
lar diameter, i.e., 



D A = - = R{t 1 )r 1 . 



(8) 



In the general case, the angular diameter distance, 
Dls{ z l, zs) = ; between two objects, for ex- 

ample, a lens at zl and a source (galaxy) at zs, reads 

H- 1 



D L s(z L ,z s ) 



(i + zs )\n k \y 2 



(9) 



dx 



x 2 f{Q, j ,x) 



III. CONSTRAINTS FROM LENSING 
STATISTICS 



is the density parameter associated with the crossover 
radius r c . For a flat universe, the normalization condition 



In this paper we work with a sample of 867 [z > 1) 
high luminosity optical quasars which include 5 lensed 
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quasars. These data are taken from optical lens surveys 
such as the HST Snapshot survey [£2], the Crampton 
survey the Yee survey p3, Surdej survey p5[ , the 
NOT Survey |§ and the FKS survey |j . Since the lens 
surveys and quasar catalogs usually use V magnitudes, 
we transform my to B-band magnitude by using B — V — 
0.2 as suggested by Bahcall et al. pqj. 

The differential probability dr of a beam having a lens- 
ing event in traversing dz^ is [§91 S3 



dr = F*(l + z L ) 2 



where 



and 
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£*ol, £*os an d Z?ls are, respectively, the angular di- 
ameter distances from the observer to the lens, from 
the observer to the source and between the lens and 
the source. For simplicity we use the Singular Isother- 
mal Model (SIS) for the lens mass distribution. The 
Schechter luminosity function is adopted and lens param- 
eters for E/SO galaxies are taken from Loveday et al. |3lf| 
(LPEM), i.e., <f>* = 3.2 ± 0.17/i 3 10" 3 Mpc" 3 , a = 0.2, 
7 = 4, = 205.3 km/s and F* = 0.010. It is worth 
mentioning that, although the recent galaxy surveys have 
increased considerably our knowledge of the galaxy lumi- 
nosity function, they do not classify the galaxies by their 
morphological type j32|. In this work we restrict our- 
selves to the LPEM parameters because they have been 
derived in a highly correlated manner and they also take 
into account the morphological distribution of the E/S0 
galaxies Q. 

The differential optical depth of lensing in traversing 
dzL with angular separation between <j> and <j> + d(j), is 
given by 
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where 4> — A6/8ir(v^/c) 2 , with the velocity dispersion u» 
corresponding to the characteristic luminosity L* in the 
Schechter luminosity function. The total optical depth 
is obtained by integrating dr along the line of sight from 
z = (zo) to 25. One obtains 

t(zs) = ^[D os (1 + z l )} 3 R 3 . (14) 
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FIG. 1: The normalized optical depth (t/F*) as a function of 
the source redshift (zs) for some selected values of Q Tc - In all 
curves, the value of the matter density parameter has been 
fixed (On = 0.3). 



Figure 1 shows the normalized optical depth as a function 
of the source redshift (zs) for Q m — 0.3 and values of 
f2 rc = 0.1, 0.2 and 0.3. For comparison, the standard 
prediction (J7 rc = 0) is also displayed. Note that, at 
higher redshifts (z > 2.5), an increase in f2 rc at fixed 
fl m tends to reduce the optical depth for lensing. For 
example, at zs = 3.0 the value of t/F* for f2 rc = 0.1 
is down from the standard value by a factor of ~ 1.10. 
This decrease of the optical depth as the value of fl Tc is 
increased (at a fixed 25 and O m ) occurs because, at high 
redshift, say z > 2.5, the distance between two redshifts 
(e.g., zo and zs) is smaller for higher values of f2 Ic . 

In order to obtain the correct lensing probability we 
have made two corrections to the optical depth, namely, 
magnification bias and selection function. Magnification 
bias, B(m, z), is considered in order to take into account 
the increase in the apparent brightness of a quasar due 
to lensing which, in turn, increases the expected number 
of lenses in flux limited sample. The bias factor for a 
quasar at redshift z with apparent magnitude m is given 

by i§ 



B(m,z) = M 2 B{m,z,M ,M 2 ), 



(15) 



where 
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In the above equation (dNq{m, z) / dm) is the measure 
of number of quasars with magnitudes in the interval 
(m,m + dm) at redshift z. Since we are modeling the 
lens by a SIS profile, Mo = 2, we adopt M 2 = 10 4 in the 
numerical computation. 

We use Kochanek's "best model" ||l[ for the quasar 
luminosity function: 



dN Q 
dm 

where 



(m, z) oc (io- a (™-" 1 ) + io-Km-w))-^ 



m 




(z - 1) for z < 1 

for 1< z < 3 
0.7(z - 3) for z > 3 



(17) 



(18) 



and we assume a = 1.07 ± 0.07, b = 0.27 ± 0.07 and 
m = 18.92 ± 0.16 at B magnitude 

The magnitude corrected probability, pt, for a given 
quasar i at Zi and apparent magnitude to, to be lensed 
is 



Pi = r(^)B(mj,^). 



(19) 



Due to selection effects the survey can only detect 
lenses with magnifications larger than a certain magni- 
tude Mf. It can be shown that the corrected lensing 
probability and image separation distribution function 
for a single source at redshift zs arc 21, |i) 

Pi(m,z) =1 
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M f = M (f +!)/(/-!) with /= io°- 4Am W. (23) 



Equation ( |21| ) defines the configuration probability, i.e., 
the probability that the lensed quasar i is lensed with 
the observed image separation. To obtain selection func- 
tion corrected probabilities, we follow f2l| and divide our 
sample into two parts, namely, the ground based surveys 
and the HST survey. 

In order to constrain the parameters £l m and Q rc we 
perform a maximum-likelihood analysis with the likeli- 
hood function given by [|TJ 



£ = Yli 1 -p'i) Y[p'kPc 



ck • 



(24) 
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FIG. 2: Confidence regions in the plane f2 m — f2 re arising from 
lensing statistics. Solid lines indicate contours of constant 
likelihood at 68% and 95.4%. 



where Nl is the number of multiple-imaged lensed 
quasars, Nu is the number of unlensed quasars, and p k 
and p % ck are the probability of quasar k to be lensed 
and the configuration probability defined, respectively, 
by Eqs. © and (||). 

Figure 2 shows contours of constant likelihood (68% 
and 95.4%) in the parameter space Q rc — Q m . The maxi- 
mum value of the likelihood function is located at Q nl = 
and fl ra — 0.03. At the la level, our analysis requires 
O ro < 0.94 and f2 rc < 0.19. Such a result means that 
a large class of these particular scenarios of brane world 
cosmology studied here are compatible with the current 
gravitational lensing data at this confidence level. In Fig. 
3a the expected number of lensed quasars, = Pi 
(the summation is over a given quasar sample), is dis- 
played as a function of Q rc with the matter density pa- 
rameter fixed at f2 m = 0.3 (as indicated by clustering esti- 
mates pql). The horizontal dashed line indicates til = 5, 
that is the number of lensed quasars in our sample. By 
this analysis, one finds tt rc ~ 0.08, a value that is very 
close to that obtained by Deffayet et at [jlTj (Q rc = 0.12) 
using SNe la and CMB data and also with the same fixed 
value for the matter density parameter. In Fig. 3b we 
show the contour for five lensed quasars in the parametric 
space fi m — tt rc . The shadowed horizontal region corre- 
sponds to the observed range fl m = 0.3 ± 0.1 We 
observe that the higher the value of fi m the higher the 
contribution of Q rc that is required to fit these data. 

At this point it is interesting to estimate the value of 
r c (the crossover distance between 4-dimensional and 5- 
dimensional gravities) from our estimates of f2 rc . In this 
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FIG. 3: a) Predicted number of lensed quasars as a function 
of f2r c for a fixed value of the matter density parameter (Q m = 
0.3) and image separation A9 < 4. b) Contour for five lensed 
quasars in the parametric space S7 m — f2 rc . The shadowed 
horizontal region corresponds to the observed range f2 m = 
0.3 ±0.1 (35). 

case, an elementary combination of our best fit (£l rc = 
0.03) with Eq. (3) provides 

r c ~ 2.8H-\ (25) 

while at the la level (il rc < 0.19), we have 

r c > 1.14H- 1 . (26) 

The former value is considerably larger than that found 
in Refs. 0, |§, i.e., r c ~ 1AH' 1 and r c ~ 0.94F" 1 



in analyses involving SNe la + CMB and angular size of 
high-z sources data, respectively. However, it is worth 
mentioning that the estimate of r c obtained in Ref. Jl7j 
refers to a flat model in which the value of the matter 
density parameter was fixed in fi ro = 0.3. As we have 
seen, by fixing this value for r2 m , the predicted number 
of lensing quasars (Fig. 3a) requires fi rc ~ 0.08 which, 
in turn, implies r c ~ 1.76H~ . 



IV. CONCLUSION 

The recent observational evidences for a presently ac- 
celerated stage of the Universe have stimulated renewed 
interest for alternative cosmologies. In general, such 
models contain an unkown negative-pressure dark com- 
ponent that explains the SNe la results and reconciles the 
inflationary flatness prediction (Six = 1) with the dynam- 
ical estimates of the quantity of matter in the Universe 
(fl m ~ 0.3 ± 0.1). In this paper we have focused our 
attention on another dark energy candidate, one arising 
from gravitational leakage into extra dimensions ^j. 
We have shown that some particular scenarios based on 
this large scale modification of gravity are in agreement 
with the current gravitational lensing data for values of 
f^in < 0.93 (la). If one fixes f2 m to be ~ 0.3, the pre- 
dicted number of lensed quasars requires Q rc ~ 0.08. 
This is a slightly open universe with a crossover radius 
of the order of r c ~ l,76H~ l . 
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